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ABSTRACT 


The  methods  of  geometrical  optics  are  extended  so  that  they  may 
be  applied  to  gyrotropic  bodies.  Various  internal  and  external  reflections 
are  considered  at  non  parallel  planar  interfaces  and  means  of  determining 
the  ray  path  or  directions  of  energy  flow  are  derived.  Non-plaiiar  geo¬ 
metries  may  be  represented  by  the  tangent  planes  at  the  various  points  of 
incidence.  A  method  is  given  for  computing  the  phases  of  the  various 
fields.  These  values  may  be  used  to  determine  the  reflected  fields  from 
such  a  gyrotropic  body. 
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GEOMETRICAL  OPTICS  FOR  GYROTROPIC  BODIES 


I.  INTRODUCTION 

A  modified  geometrical  optics  method  has  been  developed  which 
may  be  used  to  determine  the  electromagnetic  scattering  properties  of 
finite,  homogeneous,  isotropic  bodies.  The  dielectric  sphere  and 
cylinder  have  been  used  to  illustrate  the  method[l]  but  is  is  not  re¬ 
stricted  to  such  shapes.  It  has  been  postulated  thac  these  same  met¬ 
hods  may  be  applied  to  obtain  the  electromagnetic  scattering  proper¬ 
ties  for  the  gyrotropic  case.  Similar  methods  would  also  apply  to 
any  anisotropic  body.  However,  the  usual  form  of  geometrical  optics 
or  ray  methods  essential  to  the  application  of  the  modified  geometrical 
optics  method  is  restricated  to  the  isotropic  media. 

A  version  of  geometrical  optics  is  presented  which  may  be  used 
to  apply  the  modified  geometrical  optics  method  to  gyrotropic  bodies. 

It  has  been  shown  for  the  isotropic  case  that  the  reflection  and  trans¬ 
mission  coefficients  derived  using  an  infinite  planar  geometry  apply 
remarkably  well  to  curved  surfaces  with  radii  of  curvature  as  small 
as  0.  lx.  Such  coefficients  are  thus  considered  to  be  valid,  with  the 
exception  of  the  region  of  Rayleigh  scattering. 

Reflection  and  transmission  coefficients  at  a  sharp  boundary  have 
been  obtained  by  Budden[2]  and  Wait[3],  and  numerical  data  presented 
by  Yabroff[4].  Unz  has  used  a  boundary  value  solution  to  study  the  more 
general  problem  of  transmission  and  reflection  at  the  planar  boundary 
between  two  semi-finite  gyrotropic  media[5],  and  also  to  study  the 
transmission  and  reflection  due  to  an  infinite  gyrotropic  slab[6]. 

Since  there  are  generally  two  reflected  waves  and  two  transmitted 
waves  when  a  plane  wave  is  incident  from  free  space  on  an  aniso¬ 
tropic  medium  interface,  there  are  two  reflection  coefficients  and 
two  transmissions  coefficients.  These  coefficients  are  different  for 
the  two  principle  polarizations  (parallel  and  perpendicular)  as  is  true 
for  the  isotropic  case.  The  propagation  constants  in  the  gyrotropic 
medium  for  two  waves  can  be  found  from  Snells  law 

sin  6];  =  ni  sin  sin  0p^, 


where  is  the  angle  between  the  wave  normal  and  the  normal  to  the 
interface  (see  Fig.  1). 
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Consider  a  plane  wave  incident  upon  the  gyrotropic  medium  from 
free  space>  as  illustrated  in  Fig.  1.  The  two  refractive  indices*  n^ 
and  n2*  of  two  waves  in  the  gyrotropic  medium  could  be  found  from  the 
above  equation  if  and  0p  are  known;  but  the  valuer  of  0p^  and  Qp^ 
depend  on  the  orientation  of  the  static  magnetic  field  Y(Yl,Yt),  which  in 
turn  depends  on  nj  and  nz.  Solving  Maxwell's  equations  with  the  sus¬ 
ceptibility  matrix  derived  from  the  constitutive  relation  (see  Appendix 
II)  yields  the  Booker  Quartic  equation  (see  Appendix  I).  This  equation 
may  be  factored  to  find  the  component  of  the  propagation  factor  in  the 
vertical  direction  (i.  e. ,  normal  to  the  interface).  The  other  components 
of  the  propagation  factor  are  those  of  the  incident  wave  due  to  the  re¬ 
quirement  that  tangential  fields  be  matched  at  the  boundary.  Thus  the 
major  step  in  obtaining  the  direction  of  the  wave  normal  consists  of 
factoring  this  well-known  quartic  equation. 

The  four  distinct  roots  of  the  Booker  quartic  equation  are 
associated  with  the  angle  of  incidence.  Two  of  these  roots  are  associated 
with  the  sp-called  upgoing  waves  and  the  other  two  with  downgoing  waves. 
Two  upgoing  (or  downgoing)  waves  incident  from  the  gyrotropic  medium 
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on  the  "gyrotropic  rnedium-free  space”  interface  sets  up  two  downgoing 
(or  upgoing)  reflected  waves  in  the  medium  and  two  transmitted  waves 
into  free  space.  Reflection  at  a  second  interface  not  parallel  to  the 
original  one  requires  that  Booker's  equation  again  be  solved.  However, 
a  transformation  of  coordinates  yields  the  appropriate  values  of  n,  9p 
for  the  incident  wave  and  thus  reduces  the  complexity  of  the  solution. 


The  use  lon-parallel  interfaces  represents  the  appropriate 
approximation  for  determining  the  angle  of  refraction  and  the  reflection 
and  transmission  coefficients  for  two-  or  three-dimensional  bodies  with 
curved  surfaces,  such  as  the  cylinder  or  the  sphere.  The  tangent  planes 
at  the  various  points  of  incidence  provide  the  appropriate  interfaces 
which  are  used  to  represent  the  body  to  determine  these  parameters,  as 
has  been  demonstrated  for  the  isotropic  body[l  ]. 

Determination  of  the  refractive  indices,  the  reflection  and  trans¬ 
mission  coefficients,  and  the  direction  of  the  wave  normals  does  not,  how¬ 
ever,  complete  the  geometrical  optics  formulation.  It  is  also  necessary 
that  the  direction  of  energy  flow,  herein  designated  as  the  ray  path,  must 
also  be  determined,  since,  in  general,  it  does  not  coincide  with  the  wave 
normal.  Application  of  the  concept  of  stationary  phase  yields  a  relatively 
simple  expression  for  the  ray  path. 

Conservation  of  energy  now  readily  yields  the  magnitude  of  the 
fields,  with  respect  to  some  reference  point  associated  with  any  diverg¬ 
ing  ray  system  within  the  usual  restrictions  of  geometrical  optics.  The 
phase  of  these  fields  (4*)  is  readily  obtained  by  using  the  component  of 
the  vecto^propagation  constant  (rt)  in  the  direction  of  the  ray  path  (r), 
i.  e. ,  4*  ~  n  •  r.  Thus  the  fields  may  be  written  as 

where 

Aq  e^^o  is  the  field  at  the  reference  point, 

F{f )  is  the  spatial  attenuation  factor  obtained  from  the 
conservation  of  energy,  and 

I  is  the  distance  from  the  reference  point  along  the 

vector  r  {where  i  =  r  in  Fig.  4,  assuming  o  as  a 
reference  point). 

The  above  equation  will  be  discussed  in  Section  II- D. 
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IL  WAVE  mechanism  IN  A  GENERAL 
GYROTROPIC  MEDIUM 

A.  Wave  Properties  in  a  General 
Gyrotropic  Medium 

Consider  the  case  of  a  plane  electromagnetic  wave  propagating 
in  a  gyrotropic  medium  at  some  arbitrary  angle#  with  respect  to 
the  positive  z  axis.  The  usual  spherical  coordinate  system  is  used 
(see  Fig.  1)  and  eJ^^  time  convention  is  assumed.  Then 

(i)  e“jk(Six  +  SzY  +  qz)^ 

where 


k 

c 

q 

n 

r 

and 

(2)  ■  n2  =  Si^  +  S2^  +  q^ 

The  parameter  q  is  the  solution  of  a  fourth  order  algebraic  equation 
known  as  “Booker's  quarcic''[  7,8  ]  equation  and  has  four  roots.  This 
equation  is  given  in  Appendix  A.  Two  of  its  roots  (qi ,  qj )  are  char¬ 
acterized  by 

Im(qj)<  0  and  ImCqs)  <  0  . 

These  give  waves  propagating  in  the  +z  direction  and  are  known  as 
“upgoing"  waves[9].  The  other  pair,  q3,q4,  are  such  that 

Im(q3  )  >  0,  Inn(q4 )  >  0  , 


=  —  is  the  refractive  index, 

V 
_  CJ 

c 

=  velocity  of  light  in  free  space, 

=  n  cos  0p  , 

=  vector  parallel  with  axis  of  wave  normal,  and 
=  vector  from  origin  to  the  point  {x,y,  z); 
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and  propagate  in  the  -z  direction;  consequently,  they  are  known  as 
"downgoing"  waves.  This  value  of  Im(q)  indicates  the  direction  of 
energy  flow,  while  the  value  of  Re{q)  indicates  the  direction  of  the 
phase  velocity  vector  or  the  direction  of  the  wave  normal.  The 
paramaters  Si  and  S2  satisfy 


(3a) 


Si^  +  82^  =  n  sin  0 


P  ’ 


{3b) 


.  _  Si 

cos  T  -  — -  , 
n 


cos  T 


S2  . 
=  — ^  ,  and 
n 


cos  , 

n 


where  t',  and  t",  and  0p  are  the  angles  between  the  wave  normal  and 
the  respective  coordinate  axes  (x,  y,  z).  The  angle  of  the  wave  normal, 
6,  shown  in  Fig.  1  is  readily  obtained  from  Eqs.  (3a)  and  (3b)  when 
Booker's  quartic  equation  has  been  factored,  i.  e. , 


(4)  Sp  =  tan"^  ^ 

q 

In  order  to  set  up  boundary  conditions  in  as  simple  a  manner  as 
is  practicable,  it  is  necessary  to  express  all  of  the  fields  in  terms  of 
one  field  component,  in  the  manner  developed  by  Unz[5].  This  is 
done  by  first  writing  Maxwell's  equations  for  the  gyrotropic  media  as 


(5)  [S]t1oH=  -  [I]e  -  ^  P  and 

o 

(6)  [S]E=tioH, 

where 


[S]  = 

r- 

0 

q 

-q 

0 

Sz 

-Si 

‘  ^0" 

jS2 

Si 

0 

and 


[l]  is  the  identity  matrix. 
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The  constitutive  relation 


^  r 

(';  — P  =  [m]  E 

^  o 

relates  P  and  E,  where  [M]  is  the  susceptibility  matrix.  The 
matrices  [S]  and  [M]  are  given  in  Appendix  II  as  developed  by 
Budden[  lO], 

Substituting  Eqs.  (6)  and  (7)  into  Eq.  (5)  yields 


(8)  {[S]  [S]  +  [I]  +  [M]}  E  =  0. 


Equation  (8)  is  equivalent  to  three  homogeneous  equations  with  three 
field  components>  Ex»  Ey^  and  E^-  Expressing  Ex  and  Ey  in  terms  of 
Ez>  from  Eq.  (8), 


(1  -  q2  -  S2^  +  +  (Si  Sz  +  Mxy)Ey  =  -  (Siq  +  and 

(S1S2  +  +  (1  -  -  Si^  +  Myy)Ey  =  -  (Sjq  + 

Then 

(9a)  Ex  =  IT  E  =  +Hcz)(^-q^  +Myv)+(SiS2-^M^v)(S2q+M^z) 

X  z  -  {l-q2-S22+Mxx)(i-q"-Si2-fMyy)-(SiS2+M^y)(SiSi+Myx)  ^ 


and 

{9b) 


-(^-q"-Sz"+M^^)(S2q+M^,,)+(Siq+M,^)(SiS2+M^,^) 

y-'"y  (l-q2-S2HM^)(l-q2-SlHMyy)-(SiS2+M^y){SiS2+My^) 


From  Eq.  (6)  one  may  obtain 


{9c) 


H  =  !kE  E,  , 

’lo  ^ 


(9d) 


H„=3:e  =  E  ,  and 

^  ^  rin 
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(9e) 


H^=:kE  . ^ 

2-  T1  Z  T1  2 

'o  Mn 


For  each  value  of  q  (i.  e. »  q^,  q^t  q3»  q4)»  one  may  obtain  a  set  of 
values  as  Tty.,  i1y  r\^,  E^,  E^,  E^,  h^*  Hy,  and 

Unz  has  applied  equations  of  the  above  type  to  the  boundary  between 
the  gyrotropic  medium  and  free  space  to  obtain  various  reflection  and 
transmission  coefficients  for  various  interfaces [  5  ].  In  particular,  he 
has  obtained  the  reflection  and  transmission  coefficients  for  plane  wave 
incidence  from  free  space  up  to  gyrotropic  media,  from  gyrotropic 
media  up  to  free  space,  and  from  gyrotropic  media  down  to  free  space. 
Without  loss  of  generality,  the  coordinate  system  is  chosen  so  that  Si  =S  and 
$2  =  0.  This  reduces  the  complexity  of  the  solution.  These  coefficients 
are  tabulated  in  Appendix  III  for  the  convenience  of  the  reader. 

B.  Direction  of  Ray  Path 


One  method  of  obtaining  the  direction  of  the  ray  path  consists  of 
computing  the  poynting  vector  in  the  anisotropic  medium,  A  much  sim¬ 
pler  method,  given  by  Budden[l0],  consists  of  the  application  of  the 
concept  of  stationary  phase  to  compute  the  direction  in  which  the  energy 
is  traveling. 

The  phase  of  the  electromagnetic  wave  in  the  anisotropic  medium 
is  given  by 


(10)  4)  =  k(Six  +  S2y  +  qz). 


The  phase  is  stationary  when  its  first  partial  derivative  vanishes  or 


and 

(lib)  =  k  (y z)  =  0- 

9  02  9  02 
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The  ray  path  is  therefore  along  the  direction  that  satisfies  Eqs.  (H)  or 


(IZa)  tan  B 


X  0 


8si  z  8Si 


(I2b) 


0 

z  ^  ’ 


where  0g  is  the  direction  of  the  ray  path  measured  from  the  positive  z 
axis  (the  normal  to  the  interface  as  shown  in  Fig.  1). 

Without  loss  of  generality,  the  coordinate  system  is  chosen  as  before, 
so  that  Si  =  S  and  S2  =  0,  Equation  (12)  then  becomes 


(I3a)  tan  0g  =  - 

6S1  asi  Si=  S 


(l3b)  tan  0, 


_ 


’2  S2  =  0 


Booker’s  quartic  equation  is 


F(£)  =  a  q-*  +  pq^  +  Yq2  +  6q  +  e  =  0  , 


where  the  coefficients  are  functions  of  Si  and  S2 .  For  the  lossless 
case,  the  collision  frequency  is  negligible  and 


{14a)  a  =  (I-YM  +  X(Yz^-l)  =  (1-yM  +  X(£i^Y^-l)  , 


(14b)  p  2X(Si  YzY  +S2YzY  )  =  XY'(Si£i  +S2  £2  ) 

X  * 


(14c)  Y  =  -2(1-X)(C"-X)  +  2y2(c"-X)  +  X[y'-CY/+  (SiY  +S2YJ'  ] 

Z  X  "  y  ^ 

=  -2(1-X)(C2-X)  +  2y'(C^X)  +  XY'[1-C'£3'+(Si£i+S2£2)'  ]> 
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(14d)  6  =  -2tfx(SiY^Y^+S2Y^Yy)  =  -ZC%  XY"(s,<i+S2<2), 

and 

(14e)  e  =  -  C^Y^{C^-X)-  dxiSiY^i-SzYyf 

=  {l-X)iC^-xf  -  C^Y^(C^-X)  -  dXY^iSih+Szlzf  , 


where  Iz*  and  I3  are  the  direction  cosines  of  the  Y  vector  with 
respect  to  the  x,  y,  z  axes. 

Equations  (14)  are  valid  for  all  incident  angles  or  all  values  of 
S.  Thus 


(15) 


d^(q)  =  <iF(q)  _  Q  ^ 
dSi  dS2 


Differentiating  Eq.  (14)  with  respect  to  Si  yields 


(16)  i£kl=:  9Zk)  ls_+i^q^+ji  q^+ir  q^+-^a  +  i£  =  O 

dSi  8q  asi  asi  dsi  asi  asi^  asi 

There  is  a  similar  equation  for  in  which  Si  is  replaced  by  S2 


Solving  Eq.  (16)  for  and  evaluating  the  result  at  Si  ~  S  yields 

dSi 


(17) 


(18) 


tan  0 


gsi 


|£-qHlXa"  +  i^q+— 

asi^  9Si^  asi 

3F(q) 

9q 


9P  3.  av  ?  .  dd  .  ae 


tan  0 


aS; 


.  3  r  a-y  2,^-4  _ 

as2*^  '^as2  9S2 


gs. 


aF(q) 

aq 


51  =  s 

when 

52  =  0  . 
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The  derivative  of  the  coefficients  of  Booker's  equation  may  be  obtained 
from  Eqs.  {i4). 

C.  Special  Cases 

The  above  system  of  equations  gives  a  rather  complete  method  of 
determining  the  ray  path.  Certain  useful  cases  are  now  considered  to 
demonstrate  the  application  of  these  equations.  In  particular,  for 
Si  =  S,  S2  =  0,  the  derivatives  of  the  coefficients  of  Booker's  equation 
may  be  obtained  in  the  form 


(19a) 


8a 

asT 


=  0  , 


{l9b) 


f|  =  2fxl3xy^ 


1 


(I9c)  =  [4(1-X)  -  4Y^  +  2XY2(^i^  +  f3^)]  S, 

{I9d)  IjI  =  -  2£i£3XY^(1-3S^)  , 
and 

(l9e)  =  [-4(1-X)(C^-X)  +4C^Y^-  2Y^X-XY^£i^(2-4S^  )]  S, 


where  =  1  -  . 

If  the  static  magnetic  field  Y  lies  in  the  x-z  plane  {£2  =0),  then 
the  coefficients  p,  y,  5,  and  £  depend  only  on  Si.  Thus  the  partial 
derivatives  with  respect  to  S2  all  vanish  for  S2  =  0,  i.  e. , 


(20) 


aa  _  _8^  -  =; 

QSz  "  9S2  ~  aS2  "  3^  9S2 


and 
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(21) 


=  Aan  +3P  +  Zyq  +  &  =  0. 

Bq 

Several  conclusions  may  be  drawn  from  the  above  equations.  Substi¬ 
tuting  Eq.  (20)  into  (19)  yields 


tan  0 


=  0 


or 


Thus  the  ray  path  must  lie  in  the  x-y  plane  when  the  static  magnetic 
field  and  the  incident  ray  are  in  x-y  plane,  i.  e..  Si  =  S,  Sz  ~  0,  and 
iz  -  0. 

If  the  magnetic  field  is  directed  along  either  the  x  or  the  z 
axis,  then  the  four  solutions  of  Booker's  equation  are  related  by 
q3  =  -qj  and  q4  =  -qz*  These  relations  also  hold  for  normal  incidence, 
i.  e. ,  Si  =  0,  without  this  new  restriction  on  the  magnetic  field.  This 
last  case  is  of  particular  interest  since  the  ray  path  retraces  itself 
upon  reflection,  in  the  case  of  the  slab,  as  illustrated  in  Fig.  2, 


Ray  Path  And 
Wave  Normal 


Roy  Path  And 

Wave  Normol 
Z 

-  Y,=tY 


(b) 


Roy  Poth  And 
Wove  Normal 


(C  ) 


Fig.  2  -  Ray  path  and  wave  normal  of  a  wave  normally 
incident  into  an  anisotropic  medium 
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It  is  also  of  interest  to  determine  the  special  cases  for  which 
the  wave  normal  and  the  ray  path  coincide.  There  are  three  such 
cases.  If  the  static  magnetic  field  is  oriented  in  the  y  direction, 
i.  e. ,  fz  “  ^  3-nd  fi  =  =0,  then  provided  Si  =  S  and  Sz  =  0,  Eq.  (17) 

becomes 

_  9S  3S 
®gs  “  4aqS+2yq  • 


Substituting  derivatives  obtained  from  Eq.  (14)  yields 


_  [4(l-X)-4Y^]Sq^4i-4(l-X)(C^-X)44C^Y^-2Y^X]  S 

(22)  tan  0g^  =  4(l_x-Y2)q3+2[-2{l-X)(C^-X)+2Y^(C^-X)+XY2  ]q 

S 

q 

Equations  (3a)  and  (3b)  yield 

(23)  tan  0  =  -  , 

^  q 

provided  Sz  =  0*  where  0p  is  the  direction  of  the  wave  normal  with  re¬ 
spect  to  the  2  axis.  Thus  the  wave  normal  and  the  ray  path  coincide 
for  this  case  of  the  extraordinary  wave,  which  may  also  be  proved  by 
use  of  the  Poynting  vector. 

.  If  the  incident  ray  is  normal  to  the  interface,  i.  e, ,  Si  =  Sz  ~  0» 
and  the  static  magnetic  field  is  oriented  in  the  z  direction,  i,  e.,  -  i 

and  ix  and  Iz  ^  0,  then  it  can  be  readily  shown  from  Eqs,  (17)  and 
(23)  that  the  ray  path  and  the  wave  normal  again  coincide  with 
0p  =  ®gsi  "  ^  similar  result  is  obtained  when  the  static  magnetic 

field  IS  oriented  in  the  x  direction.  These  conclusions  are  summarized 
in  Fig.  3, 

D,  Electrical  Path  Length  in 
the  Anisotropic  Media 


The  electrical  path  length  is  to  be  referenced  to  the  point  of 
incidence  shown  in  Fig.  4,  which  is  also  the  origin  of  coordinates, 
i.  e. ,  X  =  y  =  z  =5  0.  Thus  the  phase  of  the  wave  at  any  point  in  the 
medium  is  given  by  Eq.  (1)  as  If  the  plane  of  incidence  is 

the  x-y  plane,  i.  e.,  Sz  =  0,  then 
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t 


^^0 »  ^0 


Mo’^0 


A'-o’  ^0 

Normal  Incidence 


P’0  *  ^0 

Norma?  Incidence 


Fig.  3  -  Three  special  cases  where  the  ray  path 
and  wave  normal  coincide 


Z 
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(24) 


6  =  k{Sx  T  qz)  =  k(n  -  r)  . 


It  is  seen  from  Fig.  4  that  this  can  be  written  in  the  form 

(25)  o  =  knr  cos  (0g  -  ^)  - 


It  should  be  noted  that  this  method  of  obtaining  the  electrical  path 
length  is  based  on  a  plane  wave  incident  upon  a  planar  interface.  This 
is  generally  true  in  the  rigorous  application  of  Snell's  law,  but  its 
application  to  other  configurations  has  proved  highly  accurate  and 
extremely  useful.  For  the  present  application  it  is  the  best  available 
approximation. 

E.  Modification  of  S  and  q  on  Reflection 
from  a  Non- Parallel  Interface 


All  of  the  equations  for  the  reflection  and  transmission  coef¬ 
ficients  at  an  interface  are  based  on  a  coordinate  system  in  which  the 
z  axis  is  normal  to  that  interface.  K  a  ray  is  transmitted  through 
this  interface  to  a  second  non-parallel  interface,  as  shown  in  Fig,  5, 
it  becomes  necessary  to  transform  our  coordinate  system  in  order 
that  previous  results  be  applicable. 

Consider  the  phase  of  Eq,  (24)  in  the  original  x,  y,  z  coordinate 
system.  The  new  coordinates  are 

(26)  x'  =  (x-a)  cos  0+  (z-b)  sin  0 
and 

(27)  z*  =  (z-b)  cos  0  -  (x-a)  sin  6 , 


where  z‘  is  normal  to  the  second  interface.  It  follows  that 


(28)  X  =  x'  cos  0  -  z'  sin  9  +  a 
and 

(29)  z  =  z'  cos  0  +  x'  sin  6  +  b  , 


1116-45 


14 


Fig,  5  -  Ray  tracing  for  non-parallel  interfaces 


and  Eq.  (24)  becomes 

(30)  <?  =  kS(x'  cos  0  -  z'  sin  0  +  a)  +  kq{|;'  cos  0  +  x'  sin  9  +  b) 
or 

(31)  4>  =  (Scos  0  +  qsin  0}  kx’  +  (q  cos  9  -  S  sin  0)  kz'  +  kSa  +  kqb. 
Recalling  that 

(32)  S  =  n  sin  0p  ;  q  =  n  cos  9p» 
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then 


(33) 

<|)  =  kn  sin(0p+0)x'  +  kn  cos(0p+9)z'  +  kn(a  sin  0p+b 

cos  Op), 

or 

(34) 

0  =  kn  sin(0p+0)x'  +  kn  cos(0p+0)z'  4  k(n  •  r)  at  O' 

f 

or 

(35) 

c|)  =  kS'x'  +  kq'z'  +  k(n  •  r)  at  O'  . 

Note  that  the  x',  z'  in  Fig.  5  for  upgoing  wave  q'  are  all  negative 
values.  In  the  following,  the  notation  Si  Sz  indicates  the  value  of 
S'  for  ray  1  and  ray  2  at  boundary  B,  respectively. 

In  the  transformed  coordinate  system 


(36)  q'  =  n  cos  +  0) 
and 

(37)  S'  =  n  sin(0p  +  0)  . 


Hence, 


(38)' 


cos(0p+6) 
cos  0p 


and 


(39) 


sin(9p-f0) 
sin  Bp 


As  the  two  interfaces  become  parallel,  0  vanishes  and  S  =  S' 
and  Si'  =  S2'  =  sin  Qj  as  before.  However,  if  0  ^  0,  then  S'  assumes 
two  different  values  for  the  two  incident  waves,  as  may  be  seen  from 
Eq.  (39)>  since  0p^  and  differ.  Recalling  that  the  coefficients 
of  Booker's  quartic  equation  are  dependent  on  S,  it  is  clear  that  two 
different  quartic  equations  must  be  factored  to  find  the  values  qi , 
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qzt  qs'i  and  in  the  primed  system  of  coordinates.  There  are  now 
two  sets  of  values  for  q'.  One  set  of  values  (ql*  'l' (4^ ) 

is  associated  with  S'l  at  point  Oi'  (shown  in  Fig,  5);  the  other  set 

‘1(3)2*  q 4 )  is  associated  with  S*2  at  point  O2.  One  of  the 
values  (q'l  in  one  set  and  q'2  in  the  other)  is  readily  obtained  from  Eq. 
(38)  for  each  set.  Thus*  the  two  sets  of  Booker’s  quartic  equations 
r  sw  may  be  reduced  tu  two  sets  of  cubic  equations,  which  must  be 
factored  to  find  the  proper  values  of  q*3  from  one  set  associated  with 
q'l  and  q'4  from  the  other  set  associated  with  q'2. 

In  other  words,  ray  tracing  is  accomplished  by  means  of  the 
equation 


(40) 


a'q'^  +  P'q'H  Y'q'2  4  £)*q*  +  e'  =  0 


for  the  prime  system,  where  b',  y\  B',  and  e'  are  functions  of  S'  and 
_.rL, — (- - ,  -  are  all  known  quantities.  For  ray  1  at  point  Oi' 

as'  as”  as'  as' 

(cee  Fig.  5),  since  qi'  is 


(41)  ai'q'3  +  (pi'  +  Qi'qi  ')q'2  +  (yi'  +  pi  'qi'  +  Ctl'  qi'")  q' 

+  {61'  +  Yl'qi'  +  Pl'qi'^  +  ai'qi'^)  =  0, 


\vh«re  Qi',  pi',  yi',  and  ^>1'  are  associated  with  Si'.  Once  the  three 
sol  tions,  q'(2)i»  l'(4)i  ’  obtained  from  Eq.  (41),  the  relation 

has  been  estalnished  that'if 


13'  <  q(4)i  > 


and  q3'  is  the  solution  we  are  seeking  for  ray  1,  on  the  assumption  that 
j  R4  and  2  R.3  are  negligible.  Once  qs'  is  obtained,  the  angle  is 
given  by 


(43) 


tan  0p,  =  i 

qi 


I 
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and  0g3  is  obtained  from  Eq.  (17). 

Similarlyj  for  ray  2  at  point  O2*  we  may  obtain  q4*  from  a  cubic 
equation  similar  to  Eq.  (41),  but  with  respect  to  Sz’  and  qz'  at  point 
O2'.  The  coefficients  of  the  Booker's  equation  may  be  replaced  by 

Ye'*  a^nd  62'  in  Eq.  (41).  The  q4'  is  chosen  by  the  relation  as 
shown  in  Eq.  (42)  (if  q2'  <  q'(i]^  .  then  q4'  >  q'(3]^).  The  technique 
applied  for  ray  2  is  the  same  as  for  ray  1.  However  q'{4)j  is  needed 
to  calculate  1*^3',  and  q’(3]^  is  needed  to  calculate  2'  R*'  • 

We  need  to  note  that  the  symbols  Si',  S2'  used  here  should  not  be 
confused  with  the  symbols  Si,  Sz  used  in  Eqs.  (17)  and  (18).  We  have 
mentioned  that  the  plane  of  the  incident  wave  is  the  x-z  plane,  hence 
we  use  S  instead  of  Si;  and  S2»  which  is  zero,  will  not  appear  in  the 
following.  When  the  incident  wave  reaches  the  boundary  A  (shown  in 
Fig.  5  ),  S  is  equal  to  sin  0  j.  Once  the  incident  wave  is  traveling  in 
the  anisotropic  medium,  it  will  split  into  two  waves.  One  of  them 
reaches  boundary  B  at  Oi'  (shown  in  Fig.  5)  as  ray  1.  The  other 
reaches  boundary  B  at  O2'  as  ray  2. 

At  boundary  B,  the  prime  coordinate  system  is  used;  hence. 


Si'  =  sin  0''j'j  (for  ray  1) 

and 

S2'  =  sin  (for  ray  2), 

where  0'ti»  are  shown  in  Fig.  5„ 

At  this  point,  the  effect  of  the  skewed  interface  has  been 
completely  determined. 


III.  SUMMARY  OF  RAY  OPTICS  FOR  ARBITRARY 
GYROTROPIC  BODIES 

All  of  the  fundamentals  needed  to  trace  a  ray  path  through  a 
finite,  homogeneous,  gyrotropic  body,  and  to  compute  the  magnitude 
and  phase  of  the  emergent  rays,  have  now  been  developed  by  con¬ 
sidering  only  planar  interfaces.  This  might  appear  to  be  a  severe 
restriction;  however,  modified  geometrical  optics  as  applied  to  a 
general  isotropic  body  is  also  based  on  the  planar  interfaces.  Snell's 
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law  and  the  various  reflection  and  transmission  coefficients  are 
all  derived  under  the  assumption  of  infinite  planar  interfaces.  Yet 
this  modified  geometrical  optics  method  has  been  remarkably 
successful  in  computing  the  radar  cross  section  for  isotropic 
spheres  and  cylinders.  This  modified  geometrical  optics  method 
also  yields  remarkably  accurate  radar  cross  sections  for  the  plasma 
cylinder  with  an  axial  static  magnetic  field[ll]. 

The  same  assumption  made  in  developing  the  modified  geometrical 
optics  method  in  these  previous  cases  would  also  be  made  here  for  this 
general  anisotropic  case.  That  assumption  is  that  the  planar  interfaces 
of  Fig.  5  would  represent  the  tangent  planes  of  any  curved  body  at  the 
point  of  incidence  and  the  point  of  reflection.  The  ray  technique  which 
has  been  developed  for  finite  homogeneous  anisotropic  bodies  will  now 
be  summarized. 

A  ray  is  incident  at  the  origin  as  shown  in  Fig.  5.  The  value  of 
S  is  obtained  from  the  incident  angle  as  S  =  sin  and  C  =  cos  0p 
where  it  is  recalled  that  the  incident  plane  wave  is  given  by 


=  Uo  e“j^(Sx+Cz)  ^ 


The  parameters  of  the  gyrotropic  plasma  medium  are  given  by 
X  and  Y  which  are  given  in  Appendix  IL  The  direction  cosines*  fi  and 
£3,  of  the  magnetic  field  vector  are  obtained  in  the  x,  y>  z  coordinate 
system.  These  represent  the  physical  parameters  of  the  anisotropic 
body  to  be  treated.  If  that  body  has  a  curved  surface  it  should  be  noted 
that  S  is  a  function  of  the  point  at  which  the  ray  enters  the  body  for  a 
particular  plane  wave  incidence;  or,  conversely,  that  the  coordinate 
system  of  Fig.  5  is  always  chosen  so  that  the  z  axis  is  normal  to  the  surface 
at  the  point  of  incidence  which  coincides  with  the  origin.  Thus  the 
values  of  S,  C,  £1,  £2,  and  £3  are  all  functions  of  this  angle  of 
incidence.  0j  . 

Once  0j  is  found,  values  of  S,  C,  £1,  £2,  and  £3  are  readily 
obtained.  These  parameters  uniquely  determine  the  values  of  the 
coefficients  of  Booker's  quartic  equation  y,  6,  e)  given  by 

Eqs.  (14a-e). 


It  should  be  recalled  that  for  our  present  case  £2  is  set  equal  to 
zero  to  maintain  the  refracted  angle  in  the  plane  of  incidence.  This  is 
not  an  essential  step  but  simply  reduces  the  complexity  of  the  solution. 
The  more  general  case  of  £2  ^0  would  follow  the  procedures  given 
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without  any  complication. 

Once  o,  y,  and  6  are  obtained,  Booker's  quartic  equation, 
Eq.  (14),  may  be  factored  ar.d  the  four  values  of  q  obtained. 

The  two  roots  qi ,  q2{for  which  the  signs  of  the  ray  path  angles 
Qgiand  9g2  the  same  as  0i)  are  the  ones  being  sought  since  they 
represent  waves  traveling  in  the  positive  z  direction.  Values  of  q^ , 
qp  may  not  always  be  positive  and  represent  waves  traveling  in 
positive  z  direction.  However,  the  ray  follows  the  path  along  which 
the  energy  flows.  This  is  a  practical  way  to  check  whether  the  qi , 

Op  are  correct. 

The  values  of  the  angles  0gj  and  Qg^  are  obtained  from  Eq.  (17). 
These  angles  are  the  directions  of  the  two  ray  paths  associated  with  qi 
and  qp,  respectively. 

The  angles  Qpj  and  of  the  refractive  index  n  will  be  useful  in 
the  following  calculations,  and  are  readily  obtained  from  Eq.  (23) 
using  values  of  qi  and  qp,  respectively.  The  values  of  m  and  n2  are 
obtained  from  Eq.  (3).  The  relative  phase  of  the  ray  at  the  point  of 
intersection  with  the  second  interface  is  obtained  by  use  of  Eq,  (25). 

It  is  necessary  to  determine  this  relative  phase  since  the  phasor  sum 
of  the  various  scattered  field  components  is  to  be  computed  to  find  the 
total  scattered  field. 

In  order  to  specify  the  magnitude  of  the  ray  at  the  point  of 
intersection  on  the  second  interface,  it  is  necessary  to  compute  the 
transmission  coefficients  at  the  first  interface.  These  are  given  by 
Eq.  (49).  At  this  time  the  reflection  coefficient,  which  would  be 
associated  with  the  direct  reflected  ray,  may  also  be  computed  at  the 
first  interface. 

The  fields  associated  with  the  two  rays  at  their  point  of  inter- 
second  interface  are  given  by 

=  Uq  "^1  ‘Sg, 


=  UoiiTz  e"j''t%r2C°s(0p2-9g2)]  , 


section  with  the 


U, 


(lA) 


and 


U 


(2A) 
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where  n  and  rz  are  the  lengths  of  the  ray  paths.  The  direct  reflected 
fields  are  given  by 


rR 


«U.f(A)  =  Uo  .,R„o 


^-jk(Sx-Cz) 


and 


rR 


-jk(Sx-Cz) 


It  is  now  necessary  to  determine  a  new  set  of  parameters  associ- 
ated  with  the  transmission  and  reflection  at  this  second  interface.  The 
first  step  is  to  transform  the  coordinate  system  such  that  the  zi'  axis 
is  now  normal  to  the  second  interface*  as  illustrated  by  the  xi',  yi'  >  zi' 
coordinate  system  for  ray  1  shown  in  Fig.  5.  The  values  qi'  and  Si'  of 
the  ray  incident  on  this  second  interface  in  this  primed  coordinate 
system  are  given  by  Eqs.  (38)  and  (39).  The  subscript  (1)  on  the 
coordinate  axes  refers  to  the  system  associated  with  qi'.  A  second 
coordinate  system  designated  by  xz’t  yzt  "z-z  is  associated  with  qz'  for 
ray  2.  This  is  necessary  because  ray  paths  ri  and  tz  differ. 

The  value  Ci'^  =1  _  Si'^  and  the  values  of  the  direction  cosines 
fi'  and  of  the  magnetic  field  are  readily  obtained  in  this  new  co¬ 
ordinate  system.  Thb  coefficients  of  Booker's  quartic  equation  (a', 

P'*  v'»  S',  £')  now  may  be  obtained  from  Eqs,  (l4a-e). 

For  ray  1*  these  coefficients  are  functions  of  Si'  at  point  Oi'; 
thus  Booker's  equation  must  be  factored  to  obtain  the  roots  qi'  q'(2)i* 

93*'  •  However  qi'  is  already  known  and  Booker's  equation 

reduces  to  a  cubic  equation:  Eq.  (41).  This  equation  is  factored  to 
obtain  the  remaining  three  roots  q*(2)i*  qs'*  q'  (4)1*  coupling  term 

is  to  be  neglected  in  the  following  work  but  there  is  no  real  need  to  do 
so  except  to  minimize  complications. 

The  desired  value  of  q'3  is  chosen  so  that  if  q'l  >  q'(2)i»  then 
93  9'(‘*)i»  and  the  associated  pair  (q'l*  9's  )  is  being  sought.  Now 

(for  ray  1)  O'gj*  as  shown  in  Fig,  5*  is  found  from  Eq,  (17 )»  and  0'pj 
from  Eq.  (43),  The  index  of  refraction,  n'3,  is  found  from  Eq,  (32), 
The  same  technique  can  be  applied  to  find  q'4  from  q'2  at  point  O2 
for  ray  2. 

The  point  of  intersection  of  these  ray  paths  with  the  next  interface 
may  then  be  found  for  the  geometry  being  treated.  In  the  example  of 
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Fig.  5j  this  is  the  original  interface.  Calculation  of  phase  at  this 
intersection  follows  the  method  described  above.  The  reflection  and 
transmission  coefficients  at  this  boundary,  i.  e.,  at  the  point  O'  of 
Fig.  5,  are  given  by  Eq.  (50).  At  point  O'l  the  values  q'l,  q'3,  q'(4)i 
are  used  to  calculate  i'R3'  and  i'T(|,  orj.)'  point  O2  the  values  q'2  » 
q'4,  q'{3)2  are  used  to  calculate  2*^4'  and  2'T(„orx)* 

The  fields  transmitted  into  free  space  at  the  second  interface 
are  given  by 


and 


"^"(iB)  i|Ti  1^11  ® 

nUljlB)  =  Uo  „Ti  e 

2  '^•1  ® 

ii^^ZB)  =  Uo,T2  2'T^e 


-Pi  g-jk{Si’xi'  +Ci'zi'), 


e-jk{S2*X2'+C2'z2')  ^ 


-j‘p2  Q-jkiSz'xz'+Cz'zz')  . 


where 


•bl  =kniri  cos(0p^-0gj) 


(pz  =  kn2r2  cos(0p2-0g2) 


and  all  phases  are  referenced  to  the  original  origin  of  coordinates  at 
point  O,  (i,  e. ,  x^  =  y  =  z  =  0),  The  fields  of  the  ray  reflected  back 
to  the  anisotropic  medium  to  the  point  Oi"  and  O2"  are 


U, 


R 


(IB) 


=  Uo  ,T,  ,R, 


and 


U 


(2B) 


=  Uo  „T 


2  2 


iRi' 


g-j(‘J>2+^2)  , 
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where 


9i‘  =  kn3  rl  cos (9^  ~  ^ 

and 

q>2  =  kn**  ri  cos(e'p^  -  0'p^  )  , 


and  all  phases  are  referenced  to  the  original  origin  O. 

By  the  use  of  ray  optics  for  the  anisotropic  bodyj  a  ray  incident  upon 
the  body  from  the  external  medium  has  been  considered  and  followed  com¬ 
pletely  through  one  internal  reflection.  All  of  the  fields  associated  with 
this  case  have  been  given;  and  any  additional  internal  reflections  may  be 
treated  simply  by  repeating  these  same  steps.  In  addition  any  coupling 
terms  may  also  be  readily  included  in  any  case  where  coupling  becomes 
significant. 

Any  changes  in  amplitude  introduced  by  diverging  ray  system s> 
i.  e.j  spatial  attenuation,  have  been  neglected.  However,  this 
problem  can  be  handled  by  the  modified  geometrical  optics  method. 
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appendix  I 

BOOKER'S  QUARTIC  EQUATION 


The  wave  form  is 


EXP{-jk(Six  +  S,y  +  qz)}  , 
and  we  may  write  symbolically 
-jkSi  , 

-  jkSi , 


_a_ 

dx 


_d_ 

dy 


and 


_3_  _ 
3z 


-jkq, 

where  hfexwell's  equation  give 
(44) 


0 

q 

'Sz 


■q 

0 

Si 


lEx 
-Si  :Ey 

Oj  E, 


H. 


‘X  ; 


0 

-q 

Sz 

iHx"! 

1 

qo 

q 

-Sz 

0 

Si 

-Si 

0 

K'  = 

Ey 

1 

"  1 

X 

Py 

_P2_ 

Combining  the  above  equations  ard  using  the  relation 

■—  Ip]  =  [m]  [eJ  , 


we  get 
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(45) 

1 

1-q  -82^4  Mxx 

S1S24  Mxy 

Siq  4  Mx2  ^ 

1 

S1S2+  Myx 

l-q^-Si^4Myy 

S2q  4  My2 

0 

II 

Siq  + 

S2q  +  M^y 

f 

! 

1. 

The  form  of  the  above 

equation  becomes 

(4  5) 

F(q)  =  Qq^  4  Pq3  + 

Yq^  4  6q  +  e  = 

0  . 

which  is  the  well-known  Booker's  quartic  equation. 

Using  the  elements  of  [m]  we  obtain 

a  =  U{U^  -  Y^)  +  X(y|  -  )  , 

P  =  2X  (SiY^Y^  +  SzY^Yy)  , 

Y  =-2U(U-X)(C^U-X)  +  2Y^(C^U-X)  + 

X  {Y^C^Yz"  +  (Si  Yx+SzYy)"}  , 

6  =  -ZC^X  (Si'^zYx+Y^Yy  )  , 

and 

e  =  (U-X)(C^U-X)^-  C^Y^(C^U-X)  -  C^X  (SiYx+SzYy)^  . 
All  symbols  X,  Y,  U  are  defined  in  Appendix  II. 
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APPENDIX  II 

THE  SUSCEPTIBILITY  MATRIX 


_  “^^N 

Ne^ 

- 

_ 

po®^o 

(0 

mw 

II 

1— ■ 

1 

N 

II 

V 

1:7 

a)=  operation  frequency, 

V  =  collision  frequency, 

cyclotron  frequency, 

=  plasma  frequency, 
m  =  mass  of  electrons,  and 


The  constitutive  relation 


-  ^o  X  E  =  UP  +  j  (Y  X  p) 


«o[E)  =  [Y]  [P]  . 


^o  ~  permittivity  in  free  space, 
p-o  =  permeability  in  free  space, 

N  =  number  of  electrons  per 
unit  volume, 

e  =  charge  of  an  electron, 

Hq  =  applied  static  magnetic 
field. 


where 
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(47) 


gives 


U 


-JYz  jYy  ] 

U  -jYx  I 


[-JYy  jY^  U 


I 


-L  [P]  =  [Y]-‘  [E]  =  [M]  [E]  , 


where 


(48) 


[M] 


-X 

U(U^-Y^) 


U"-Yx  jUY^-Y^Y. 

-jUYz-YxYy  U^-Yy^ 


jUYy-YxY^  -jUYx-YyY 


and  P  is  electric  polarization. 


-jUYy-Y^Y 

jUY^-YyY 

U^-Yz^ 
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APPENDIX  III 

REFLECTION  AND  TRANSMISSION  COEFFICIENTS 
AT  PLANAR  "FREE-SPACE  GYROTROPIC"  INTERFACES 


The  coordinate  system  used  in  this  Appendix  is  chosen  so  that 
Si  =  S  and  Sj  =  0.  There  are  four  possible  waves  in  the  gyroptropic 

medium;  paired  as  upgoing  waves,  designated  by  subscripts 
(1)  and  (2),  and  downgoing  waves,  designated  by  subscripts  (3)  and  (4), 
Thus  a  wave  (1)  or  (2)  in  the  gyrotropic  medium  may  be  reflected  as 
wave  (3)  or  (4),  the  reflection  coefficient  of  which  is  designated  i  R3  , 
iR4,2F3»  or  aUt*  etc.  The  quantities  r\,  ir  are  obtained  from  Eq.  (9) 
by  setting  S2  =  0. 

A.  Wave  Incident  From  Free  Space  up 
to  Gyrotropic  Medium 

For  parallel  polarized  plane  wave  incidence  (Ey  =  0) 


T2  = 


E 

=  C — b- 


Di  c  ^xi)  » 


"  ffir  ~  ■%~r  ~  •^t^^XlV2"%l^X2)  ■^•Q^^Xl^X2"'nxi'"‘x2) 
y  ^ 

C('iryjriy2 -TlyjTTy^  )  -  (^Xl’*^y2  ”^yi^X2  ^ 


""  IoHJ  ^  1 


For  perpendicular  polarized  plane  wave  incidence  (Ex^~ 


(49b) 
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{49b) 

(cont) 


X  Tz  =  -  £)j  (WXI  +  C^yj  )  » 


TT  R  1 

x^x  ^  =-^  f  "(’Ixi^yz  “^yi’^xa^  ■*■  ~  ^’^xi^xz  "^^xi^xz  ^ 

Ey^  Di  C 

C(Tryjiiy.^  “’^yi^yz  ^  ^  ^’^xi’^^yz  ”  ’^yi^^xz  ^ 

and 

Tl  H  ^  Z 

X^l  “  ^  (’^XlVz  "\l^X2  ^  » 

where 


1 

"  ^’^xi'ny2  “^yi’^X2  ^  “  (3  ^^xi^X2  “^xi^X2  ^ 


C(Tryjiiy^-TiyiTry^)  +  (^xi^y2"Vl’*^X2^  * 


Transmission  into  the  gyrotropic  medium  yields  coefficients  of  the 
form  |,Ti ,  II  T2  »  v/here  the  subscript  it  designates  the  polarization  of  the 
incident  wave  and  (l)or  (2)  designates  one  of  the  upgoing  waves.  Conven¬ 
tional  definitions  of  polarization  are  used  so  that  subscript  11  or  x  means 
that  the  E  vector  is  "parallel  to"  or  is  "perpendicular  to"  the  plane  of 
incidence*  respectively. 


B.  Wave  Incident  From  Gyrotropic  Medium 
Up  to  Free  Space 


For  wave  (1)  (E2;2^  =  0) 


(50a) 


^"^11  ~  E  f 


E  T 

y~  =  -1  X 


tiHJ  1 


C  E  I 
zi 


_£ 

_C 

Dz 


{^Xl(%3V-V3^X4)  -  ^Xl(Trx3V^yS»X4) 


4  ^yi(^X3^X4“^X3^X4^  ^’^xi{^y3^y4“^y3Ty4  ) 


-  CTryi(irx3Tly^-Tly3Wx4)  +  CTlyiCTrxjTTy^-WyjTfx^  )}, 
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t 


(50a) 

(cont) 


E  1 

-  ^Xl(^y3V4-’1y3’ry4) 

+  V^’^yS^X4-^X3V)  -  ~^Xl(^y3^X4-%V4^ 

+  -i  1 

C  ^yi^^X3^X4-’1x3Tx4J  "  C^Xl(^X3^y4-^y3^X4)}* 


_  R  1  ^ 

E^jj!  dJ  ^  ^^^1%-*” ^71^x4)  ■^(^xi^X4“^xi^X4  ^ 

+  C(TryjTiy^-Tiyj^y^)  +  (TTx^Try^-TyjTrx^)]  , 

and 

E2jI  ^^^xi^y3~%l^X3^  ■^^^xi^X3“^xiTrx3) 


^  ^^Vl^y3~\l^y3^  ^  ^’^Xl^y3~’‘’yi^X3)]* 


For  wave  (2) 


(50b)  2T„ 


(Ezil  =  0) 

=  j  i 

C  Dl 


^^X2(^X3^y4-^y3’lx4)  -  X^’^X3^y4-^y3V) 


+  %2  (’^X3nX4-^X3^X4)+ 

-  C  TTy^  {’rx3^y4-^y3’rx4)  +  C  ’ly2K3^y4-ffy37rx^)}, 


2-^11 


.  V 


E^l  Dz  ^^^2(^X3^y4-'ny3^X4)  -  "^Xz  ^’^y3’^y4"'^y3^y4) 


'*'  \2  ^^y3^X4~’^X3^y4)  "  C ’^Xa  ^  V3^X4-'nX3^y4) 

C  V2^^X3  V^X3%4)  -  C’^X2(’^X3V4"’^y3^X4))  » 
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(50b)  2R3  f  ^’^X2\4~\2^X4^  C  ^’x2’1x4"’^X2’^X4^ 

(cont)  ^ 

^  ^’^X2‘^y4"’’^y2’’’x4^3  * 

and 


I'Z  ^1  1 

2R^  =  ^  [ K^ly^-lyaV^)  +  -Ft’^x^lx,- Vx,) 


'Zz 


‘y2’^y3^ 


V" 


^zVi  yzX2‘ 


where 


O2  =  (Ixjlyr’lyj’lx,)  +  C  (’■xj’tx.-’lxjX^^)  +  C(X^ny,-’ly3Ty,) 


■*■  {Xx3^y4"*y22rx4) 


The  dominant  reflection  coefficients  are  1R3  and  2R4»  Coefficients  1R4 
and  2R3  are  usually  negligible. 

C.  Wave  Incident  From  Gyrotropic  Medium 
Down  to  Free  Space 

For  wave  (3)  =  0) 

•ppj  T  E  ~C 

(51a)  3  T„  =  J-  =  -  C  J  =  {^X3{^Xl'ny2-'nyiTlx2  ^  “  ^X3(’^Xl^y2"’^yi^X2) 

■^23  ^ 

C 

+  ^y3Kl^X2"^Xl’^X2)  -  C  ^X3(V^y2-^yiV^ 

+  CTry3(irxi  r!y2  -riyj^Trx2)  -  C^y3(^xi^y2  “■"■yi'n-x^)}! 
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T 


(Sla) 

(cont) 


EyT  1 

3 "El-  ”  -  pj  ^’^y3(’^xi^y2~’lyi^X2^  "  ^X3(Vl^y2"^yi^y2) 

+  \3<^i’’x2-’’x,\2’  ■5’'x3('yi’'x3->lx,’ryP 

1  ,  i 


E  ^ 

^  =  ^f(''x3\3-’V^’'x3)  -  -i  ixx3’lx3-’lx3’'x3) 
CCTTy^Tly^-tl  TTyj)  +  {^X2’’’y3"^y2^X3)} 


Ez  R  2  j^ 

"  "Efr  =  ^  {-Ki%3- V^Xj)  +  c  (^X1^X3-^X1^X3) 


+  c  (^yi’1y3-\jTTy^)  -  {Trxi'>ryj-TryjTrx3)}< 


For  wave  (4)  (Ez,^  =  0) 


(51b) 


_^Hy 


E„  I 

24 


E24I 


■  Ds  ^’'x^f’lxi'lyj-'lyi'lxj) 


•\,K\2-\2’'x3>  +  '’y4(’'xi’lx3-\/x3) 


-  <^”x4(’'y,''y3-’’y.’'y3>  +  CXy^fx^.ny^-'iy^ 


-  C\,(irx,Xy^-XyjX^^ 


■  D3  ^''y4(''xi''y2-''yi’'x2’  ■  Wv'*y2‘'’yi’'y2' 


+  ’ly4(’'y,’lx2-1xiXy2)  +  C’'x4(’'y,’lx2-’lxi’'y2> 
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{51b) 

(cent) 


-  Jl  ^  trr  W 

C  ^  C  ^X4^Xl^y2"^yi’^X2'  ■’  * 


and 


E  ^  1  1 

•<Rl  ^  ^f’'x2''yi'’ly2’'x4)  ■  ■q('x2’'x4"''x2’'x4) 


-  C  (”y2''y4-''yz’'y-.’  (’'X2’fy4-’'y2’'x4)}  ' 


E  ^  1  1 

<^2  =  EV=D:{-(’'xi''y4-’lyl1x4)  +-ct’'xi''X4-*lxi’'x4) 

Z4  ^ 


^ •^yi’iyr’'y> V)  -  (’rxiV’^yi’'^'*’^’ 


where 


^  =  (^Xl^y2-^yi%2^  -  ■^Ki^X2-^Xi’»-X2)  -  CCTTyjTly^-TlyjTry^) 

+  (’'^xi^y2“^yi’'’^X2^  * 

In  this  case  the  dominant  reflection  coefficients  are  3R1  and  4R2.  The 
coefficients  3R2»  4R1  are  usually  negligible. 
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